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■ Abstract 

We study integrability of the Euler-Poisson equations describing the motion of a rigid body 
with one fixed point in a constant gravity field. Using the Morales-Ramis theory and tools 
of differential algebra we prove that a symmetric heavy top is integrable only in the classical 
1 cases of Euler, Lagrange, and Kovalevskaya and is partially integrable only in the Goryachev- 

Chaplygin case. Our proof is alternative to that given by Ziglin (Funktsional. Anal, i Prilozhen., 
17(l):8-23, 1983; Funktsional. Anal, i Prilozhen., 31(1):3-11, 95, 1997). 

> : 

J^j ■ Resume 

Nous etudions l'integrabilite des equations de Euler-Poisson qui decrivent le mouvement d'un 
solide rigide avec un point fixe dans un champ gravitationnel constant. En utilisant la theorie 
de Morales-Ramis et des outils d'algebre differentielle, nous prouvons qu'un solide symetrique 
est integrable seulement dans les cas classiques d'Euler, Lagrange et Kowalevski, et est par- 
tiellement integrable seulement dans le cas Goryatchev-Tchaplygin. Notre preuve est une al- 
ternative a celle donnee par Ziglin (Funktsional. Anal, i Prilozhen., 17(l):8-23, 1983; Funktsional. 
Anal, i Prilozhen., 31(1):3-11, 95, 1997). 
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1 Equations of motion and motivation 



Equations of motion of a rigid body in external fields are usually written in a body 
fixed frame. Here, we use the following convention. For a vector v we denote by 
V = [Vi, V2, V^] T its coordinates in a body fixed frame, and we consider it as a one 
column matrix. The vector and scalar products of two vectors v and w expressed in 
terms of the body fixed coordinates are denoted by [V, W] and (V, W), respectively. 

We consider a rigid body of mass m located in a constant gravity field of intensity 
g. One point of the body is fixed. The distance between the fixed point and the mass 
centre of the body is D. Assuming that gD 7^ 0, we choose units in such a way that 
\x := mgD = 1. The Euler-Poisson equations 

— M= [M,JM] + [N,L], — N=[N,JM], (1) 

describe the rotational motion of the body. In the above equations M denotes the an- 
gular momentum of the body, N is the unit in the direction of the gravity field, and L 
is the unit vector from the centre of mass of the body to the fixed point; J is the inverse 
of the matrix of inertia, so O := JM is the angular velocity of the body. The principal 
moments of inertia A, B, and C are eigenvalues of J" 1 . For our further consideration it 
is important to notice that in (|TJ) the body fixed frame is unspecified, so we can choose 
it according to our needs. A body fixed frame in which J is diagonal is called the princi- 
pal axes frame. This frame is uniquely defined (up to the numbering of the axes) when 
J has no multiple eigenvalues. If J has a multiple eigenvalue, e.g. when A = B, we say 
that the body is symmetric. Then the principal axes frame is defined up to a rotation 
around the symmetry axis. 

System (|TJ) depends on parameters A, B, C, and L but physical constraints restrict 
the allowable values of parameters to a set V C M 6 defined by the following conditions 

A>0, B>0, C>0, (L,L) = 1, 
A + B>C, B + C > A, C + A > B. 

Euler-Poisson equations possess three first integrals 

H=1(M,JM) + (N,L), (2) 
Hi = (M,N), H 2 = (N,N). (3) 

It is known that on the level 

M x = {(M,N) G M 6 I Hi = X ,H 2 = 1} 

the Euler-Poisson equations are the Hamiltonian ones, see e.g. 

Remark 1. Configuration space of a rigid body with a fixed point is the Lie group 
SO (3, M) (all possible orientations of the body with respect to an inertial frame). Thus, 
classically the phase space for the problem is T*SO(3, R). Hence it is a Hamiltonian 
system with three degrees of freedom and it possesses one additional first integral Hi 
(first integral H 2 is identically equal to one in this formulation). The existence of this 
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first integral is related to the symmetry of the problem (rotations around the direction 
of the gravity field) and this allows to reduce the system by one degree of freedom. 
The Euler-Poisson equations can be viewed as an effect of reduction of the system on 
T*SO(3, R) with respect to this symmetry The phase space of the reduced problem 
can be considered as the dual g* to Lie algebra g of group of rigid motions G = R 3 x 
SO (3,1ft). We can identify (N, M) as element of q* = R 3 * x so(3,R)* using standard 
isomorphisms between so (3, R) and R 3 , and between R 3 * and R 3 . Let (X, A) G G, 
with X <E R 3 and A <E SO (3, R). Then the coadjoint action Ad* x A) : g* -> g* is defined 
by 

Ad (x,A)( N ' M ) = ( AN ' [ X ' A N] + AM). 
Functions Hi and Hi given by ® are orbits invariant, i.e., on each orbit of coadjoint 
action 

O iNM )-= U {(AN,[X,AN]+AM)} 

(X,A)eG 

they have constant values. They are called Casimirs. Orbits of coadjoint action defined 
above coincide with M x which is diffeomorphic to T§ 2 . As it is well known, orbits 
of a coadjoint action are symplectic manifolds equipped with the standard Kostant- 
Berezin-Souriau-Kirillov symplectic structure [18; 44]. On four dimensional orbits M x 
the Euler-Poisson equations are Hamiltonian with H given by (|2j) as the Hamilton func- 
tion. 

Thus, as a Hamiltonian system on M x the Euler-Poisson equations have two de- 
grees of freedom and are integrable on M x if there exists an additional first integral H3 
which is functionally independent with H on M x . Equivalently we say that the Euler- 
Poisson equations are integrable if there exists a first integral H3 which is functionally 
independent together with H, Hi and Hi. 

We say that the Euler-Poisson equations are partially integrable if they are inte- 
grable on Mo- 

The known integrable cases are the following: 

1. The Euler case (1758) corresponds to the situation when there is no gravity (i.e. 
when \i = 0) or L = (the fixed point of the body is the centre of mass). The 
additional first integral in this case is the total angular momentum H3 = (M, M) . 

2. In the Lagrange case [26] the body is symmetric (i.e. two of its principal moments 
of inertial are equal) and the fixed point lies on the symmetry axis. The additional 
first integral in this case is the projection of the angular momentum onto the 
symmetry axis. If we assume that A = B, then in the Lagrange case Li = Li = 0, 
and H 3 = M 3 . 



3. In the Kovalevskaya case |2£jt |21Q the body is symmetric and the principal mo- 
ment of inertia along the symmetry axis is half of the principal moment of in- 
ertia with respect to an axis perpendicular to the symmetry axis. Moreover, the 
fixed point lies in the principal plane perpendicular to the symmetry axis. If 
A = B = 2C, then (after an appropriate rotation around the symmetry axis) we 
have in the Kovalevskaya case Li = L3 = 0. The additional first integral has the 
form 

H 3 = Q(M? - Ml) + Ni\ + (M1M2 + N 2 ) 2 . 
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4. In the Goryachev-Chaplygin case 115] the body is symmetric and, as in the Ko- 
valevskaya case, the fixed point lies in the principal plane perpendicular to the 
symmetry axis. If we assume that the third principal axis is the symmetry axis, 
then in the Goryachev-Chaplygin case we have A = B = 4C and L2 = L3 = 0. 
In the Goryachev-Chaplygin case equations ((TJ) are integrable only on the level 
Hi = and the additional first integral has the following form: 

H 3 = M 3 (M? + M\) + M1N3. 



For a long time the question if there are other integrable cases of the Euler-Poisson 
equations except those enumerated above was open, although many leading mathe- 
maticians tried to give a positive answer to it. 

The problem was completely solved by S.L. Ziglin in a brilliant way. First, in |0] 
he proved the following. 

Theorem 1 (Ziglin, 1980). If (A — B)(B — C)(C — A) 7^ 0, then the Euler-Poisson system 
does not admit a real meromorphic first integral which is functionally independent together 
with H, Hi and Hi. 



Later, he developed in [50] an elegant method, which now is called the Ziglin the- 
ory, and using it he proved in |5l|l the following theorem. 

Theorem 2 (Ziglin, 1983). The complexified Euler-Poisson system for a symmetric body is 
integrable on Mo with complex meromorphic first integrals only in the four classical cases. 

This result was improved in [52]. 

Theorem 3 (Ziglin, 1997). The Euler-Poisson system for a symmetric body is integrable on 
Mq with real meromorphic first integrals only in the four classical cases. 

We know that in the Euler, Lagrange and Kovalevskaya cases the Euler-Poisson sys- 
tem is integrable globally, i.e. the additional first integral exists on an arbitrary sym- 
plectic manifold M x . In the Goryachev-Chaplygin case the Euler-Poisson equations 
are integrable only on Mo, and this fact was also proved by Ziglin. 

Remark 2. We do not even try to sketch the very rich history of investigations of the 
problem of the heavy top. We refer here to books [14; 27; 23; S] and references therein. 
However, several works where the question of integr ability of the problem was inves- 
tigated are worth mentioning. In [22 ] it was shown that for a non-symmetric body the 
Euler-Poisson equations do not admit an additional real analytic first integral which 
depends analytically on a small parameter \i, see also [23, Ch. III]. For a symmetric 
body when the ratio of the principal moments is small enough (i.e. for the case of the 
perturbed spherica l pe ndulum ) the non-existence of an additional real analytic first in- 
tegral was proved $24; 25 ]. A similar result for a perturbed Lagrange case was shown 
in fli: 0/. A novel, variational approach to proving the non-integrability was elabo- 
rated by S. V. Bolotin in £3/ where he showed the non-existence of an additional real 
analytic first integral for a symmetric heavy top for the case when the fixed point lies 
in the equatorial and the ratio of the principal moments of inertia is greater than 4. 
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The Ziglin theory is a continuation of the idea of S. N. Kovalevskaya who related 
the (non)integrability with the behaviour of solutions of the investigated system as 
functions of the complex time. The main object in the Ziglin theory is the monodromy 
group of variational equations around a particular non-equilibrium solution. As it was 
shown by Ziglin, if the investigated system possesses a meromorphic first integral, the 
monodromy group of variational equations possesses a rational invariant. Thus, for 
an integrable system, the monodromy group cannot be too 'rich'. The main difficulty 
in application of the Ziglin theory is connected with the fact that 'except for a few differ- 
ential equations, e.g., the Riemann equations, Jordan-Pochhammer equations and generalised 
hypet -geometric equations, the monodromy group has not been determined' , [45, p. 85]. For 
other second order differential equations only partial results are known, see e.g. 0:0]. 
Having this in mind, one can notice that the analysis of the variational equations and 
determination of properties of their monodromy group given by Ziglin in his proof of 
Theorem 12 is a masterpiece. Later, the Ziglin theory was developed and applied for 
a study of non-integrability of various systems but, as far as we know, nobody used 
Ziglin's brilliant technique developed in his proof of Theorem |2 

In the nineties the theory of Ziglin was extended by a differential Galois approach. 
It was done independently by C. Simo, J. J. Morales-Ruiz, J.-P. Ramis 1E7L I38LI39L Eoll 
and A. Braider, R. C. Churchill, D. L. Rod and M. F. Singer O Q. Nowadays, this 
approach is called the Morales-Ramis theory. The key point in this theory is to re- 
place an investigation of the monodromy group of variational equations by a study 
of their differential Galois group. The main fact from this theory is similar to that of 
Ziglin: the existence of a meromorphic first integral implies the existence of a rational 
invariant of the differential Galois group of variational equations. Forgetting about 
differences in hypotheses in main theorems of both theories, the biggest advantage of 
the Morales-Ramis theory is connected with the fact that applying it, we have at our 
disposal developed tools and algorithms of differential algebra. Thanks to this fact, it 
can be applied more easily. 

We app lied the Morales-Ramis theory to study integrability of several systems, see 
e.g. 112 81: l29c l32c l30r 13 it l34r 13311 , and we notice that obtaining similar results when work- 
ing only with the monodromy group is questionable, or, at least difficult. This gives 
us an idea to reanalyse the Ziglin proof of Theorem |2 which is rather long (about 10 
pages in l5llB . We wanted to present a new, much shorter and simpler proof which is 
based on the Morales-Ramis theory and tools from differential algebra. In fact, at the 
beginning, we believed that giving such proof would be a nice and simple exercise but 
quickly it appeared that we were wrong. A 'naive' application of the Morales-Ramis 
theory leads quickly to very tedious calculations or unsolvable complications. As our 
aim was to give an 'elementary' proof, we put a constrain on the arguments which are 
allowable in it: no computer algebra. Thus we spent a lot of time analysing sources of 
difficulties and complications, and the aim of this paper is to present our own version 
of proofs of Theorem |2 and Theorem |3l As we believe, these proofs present the whole 
power and beauty of the Morales-Ramis theory. 

The plan of this paper is following. To make it self-contained in the next section we 
shortly describe basic facts from the Morales-Ramis and Ziglin theory. We collected 
more specific results about special linear differential equations in Appendix. In Sec- 
tion |3] we derive the normal variational equation. Sections HI and |5] contain our proofs 
of Theorem |2 and IS In the last section we give several remarks and comments. 
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2 Theory 



Below we only mention basic notions and facts concerning the Ziglin and Morales- 
Ramis theory following |50tl4c 38]. 

Let us consider a system of differential equations 

— x = v(x), t G C, x G M, (4) 

defined on a complex n-dimensional manifold M. If <p(f ) is a non-equilibrium solution 
of (HJl, then the maximal analytic continuation of <p(t) defines a Riemann surface T 
with t as a local coordinate. Together with system (|4|) we can also consider variational 
equations (VEs) restricted to TpM, i.e. 

I = T(v)L l g T r M. 

We can always reduce the order of this system by one considering the induced system 
on the normal bundle N := TrM /TV of V [50] 

f[ = 7^(T(v)7i- l Vi), 77 EN. (5) 

Here n : TpM — ► N is the projection. The system of s = n — 1 equations obtained in this 
way yields the so-called normal variational equations (NVEs). The monodromy group 
Ai of system © is the image of the fundamental group 774 (P, to) of V obtained in the 
process of continuation of local solutions of © defined in a neighbourhood of to along 
closed paths with the base point to. By definition, it is obvious that Ai C GL(s, C). A 
non-constant rational function f(z) of s variables z = [z\, . . . , z s ) is called an integral 
(or invariant) of the monodromy group if f(g ■ z) = f(z) for allg G M. 

From the Ziglin theory we need the basic lemma formulated in [50] and then given 
in an improved form in 15211 . 

Lemma 1. If system ® possesses a meromorphic first integral defined in a neighbourhood 
U C M, such that the fundamental group of V is generated by loops lying in U, then the 
monodromy group M. of the normal variational equations has a rational first integral. 

If system (HJ) is Hamiltonian, then necessarily n = 2m and M is a symplectic man- 
ifold equipped with a symplectic form w. The right hand sides v = of © are 
generated by a single function H called the Hamiltonian of the system. For given H 
vector field Vh is defined by cv(vh, u) = dH ■ u, where u is an arbitrary vector field on 
M. Then, of course, H is a first integral of the system. For a given particular solution 
<p{t) we fix the energy level E = H(<p(t)). Restricting ® to this level, we obtain a well 
defined system on an (n — 1) dimensional manifold with a known particular solution 
(p(t). For this restricted system we perform the reduction of order of variational equa- 
tions. Thus, the normal variational equations for a Hamiltonian system with m degrees 
of freedom have dimension s = 2(m — 1) and their monodromy group is a subgroup 
of Sp(s, C). 

In the Morales-Ramis theory the differential Galois group Q of normal variational 
equations plays the fundamental role, see [38; 39]. For a precise definition of the dif- 
ferential Galois group and general facts from differential algebra see E EH USB 
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We can consider Q as a subgroup of GL(s, C) which acts on fundamental solutions of 
^ and does not change polynomial relations among them. In particular, this group 
maps one fundamental solution to other fundamental solutions. Moreover, it can be 
shown that M. C Q and Q is an algebraic subgroup of GL(s, C). Thus, it is a union 
of disjoint connected components. One of them containing the identity is called the 
identity component of Q and is denoted by Q°. 

Morales-Ruiz and Ramis formulated a new criterion of non-integrability for Hamil- 
tonian systems in terms of the properties of G° ll38l:l39il . 

Theorem 4. Assume that a Hamiltonian system is meromorphically integrable in the Liouville 
sense in a neigbourhood of the analytic curve V. Then the identity component of the differential 
Galois group ofNVEs associated with V is Abelian. 

In most applications the Riemann surface V associated with the particular solution 
is open. There are many reasons why it is better to work with compact Riemann sur- 
faces. Because of this, it is customary to compactify V adding to it a finite number 
of points at infinity. Doing this we need a refined version of Theorem HJ for details 
see [38; 39]. However, in the context of this paper, the thesis of the above theorem 
remains unchanged if instead of V and the variational equations over V we work with 
its compactification. 



3 Particular solutions and variational equations 

To apply the Ziglin or the Morales-Ramis theory we have to know a non-equilibrium 
solution. Let us assume that the fixed point is located in a principal plane. Then, in fact, 
we can find a one parameter family of particular solutions which describe a pendulum 
like motion of the body. For a symmetric body this assumption is not restrictive (if 
necessary we can rotate the principal axes around the symmetry axis). 

We choose the body fixed frame in the following way. Its first two axes lie in the 
principal plane where the fixed point is located and the first axis has direction from the 
fixed point to the centre of mass of the body. We call this frame special. A map given 

(M, N) = (RM, RN), R e SO(3, R), 
transforms equations © to the form 

= [MJM] + [N,L], = [N,JM], 



where 



J = R T JR, L = R T L. 



Now, if symbols with tilde correspond to the principal axes frame, then, taking into 
account our assumption about the location of the fixed point, we have 



A" 1 





o " 







zr 1 





L 








C" 1 





[Li,0,L 3 ] 
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Taking 



R 



Li L 3 
1 J) 
-L 3 Li. 



we obtain 



where 



[1,0, 0] T , J 



T 2 L 2 
AC' 

2^= '^"x) LlL3 ' 



' fl 2rf' 
b 
2d c 



7 2 L 2 

Zi _l Z3 

C A ' 



»4 



(6) 



Thus, the prescribed choice of R corresponds to the transformation from the special 
to the principal frame. Without loss of generality we can put b = 1. For a symmetric 
body we assume that A = 8 / C. Under this assumption, if d = 0, then I4L3 = 0, 
and, in this case, the special frame coincides with the principal frame. 

From now on we consider the complexified Euler-Poisson system, i.e. we assume 
that (M, N) G C 6 . 



3.1 Case d ^ 

It is easy to check that manifold 

M := {(M,N) G C 6 I Mi = M 3 = N 2 = 0,N 2 + = 1} C M , 

is symplectic sub-manifold of Mo diffeomorphic to T §^ C T S 2 . (by §^ we denote 
m-dimensional complex sphere). Moreover, J\f is invariant with respect to the flow 
generated by (|TJ|. The Euler-Poisson equations restricted to J\f have the following form 

A A A 

— M 2 = N 3 , — N! = -M 2 N 3 , —N 3 = M 2 N 1 , (7) 

and are Hamiltonian with H ^ as the Hamiltonian function. For each level of Hamil- 
tonian = e := 2fc 2 — 1, we obtain a phase curve 1\. We restrict our attention to 
curves corresponding to e G (—1, 1] so that G (0, 1]. A solution of system Q lying on 
the level H\ M = 2k 2 - 1 we denote by (M 2 (f, k), Ni (t, k), N 3 (t, k) ). 

For a generic value of k phase curve r\ is an algebraic curve in C 3 {M 2 ,Ni,N 3 },and, 
as intersection of two quadrics 

2k 2 - 1 = ^M 2 + Ni, N 2 + N 2 = 1, 

is an elliptic curve (for k = 1 it is a rational curve). We can compactify it adding two 
points at infinity which lie in directions (0, ±i, 1). Thus, a generic r\ can be considered 
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m\ 




2dM 2 





(c - 1)M 2 













m\ 


m 2 



















1 




m 2 


m 3 




(1 - a)M 2 





-2dM 2 





-1 







m 3 


«i 







-N 3 











-M 2 




ni 


n 2 




aN 3 - 2dN\ 





2dN 3 - cNi 













n 2 


K3 







Ni 





M 2 










n 3 



as a torus with two points removed. In our further consideration we work with sub- 
family r\ with k G (0, 1). Only in Section|5]we refer to the phase curve Pi corresponding 
to Jfc = 1. 

Equations (0) describe the pendulum-like motions of the body: the symmetry axis 
of the body remains permanently in one plane and oscillates or rotates in it around the 
fixed point. 

For a point p = (M, N) £ Mo by v = (m, n) we denote a vector in T p Mq. Varia- 
tional equations along phase curve 1\ have the following form 



d_ 
df 



where (M 2/ Ni / N 3 ) = (M 2 (t / k),N 1 (t / k) / N 3 (t / k)) G 1\. They have the following first 
integrals 

h = M 2 m 2 + ni, hi = N\tn\ + N 3 m 3 + M 2 n 2/ h 2 = N\n\ + N 3 n 3 . 
As it was shown by Ziglin, the normal variational equations are given by 

= 2dM 2 mi + (c - l)M 2 m 3 , 

~^j m 3 = (1 — a)M 2 m\ — 2dM 2 m 3 — n 2 , 

^jn 2 = (aN 3 - 2dNi)nii + (2dN 3 - cNi)m 3/ 
= Nirai + N 3 m 3 + M 2 n 2 . 

We assume that the particular solution is not a stationary point (M 2 (t, k) = 0, Ni(t, k) = 
±1, N 3 (t, k) = 0). Under this assumption we reduce the above system to the form 



df 
d_ 

dt 



m\ = 2dM 2 m\ + (c - l)M 2 m 3 , 



m 3 



Ni 
W 2 



+ (1 - a)M 2 



mi + 



N3 
M 2 



2dM 



m 3 . 



We can write the above system as one second order equation 

d 2 , . d , . 

— -~m + a\(t)— m + anitjm = 0, m = m\, 

at 2 - at 



with coefficients 



ai(t) 



N 3 (t,k) 
'M 2 {t,k)' 



(8) 



(9) 



a (t) 



c)N 1 {t,k)+2dN 3 (t,k) +fM 2 (t,k) 2 , 



and 



f =(a- l)(c- 1) -4d 2 . 
9 



Now, we make the following transformation of independent variable which is rational 
parametrisation of the complex circle §^ 



-> z :- 



N 3 (t,k) 



Then we obtain 



l + z 2 ' 



z 2 



N 3 : 
1 



2z 



l + z 2 ' 



M 2 

„2\ 



2z 



l + z 2 



where 



1 



1+e 



I' 



After transformation ((TO)) equation © reads 

ra" + p(z)tn' + q{z)m = 0, 

with coefficients 



p(z) 



- + 



1 



_d 
dz' 

1 



z + i 



z + s 



and 



<?(z) 



4 

E 



OCi 



where we denote [z\, Zj_, Z3, Z4} 



& -i, s, 



«i 



(5 



c) - f + id, 
2d 



ft 



1 + s 2 
l-c)(l-s 2 )+4ds 
2s(l + s 2 ) 



+ 

-s), and 

0C2 = oc\, 
c-1 

TTs 2 

/3 4 = 



ft 



0C3 ■ 
Pi 



«4 = 0, 



01, 



+ Ads 



2s(i+s 2 : 



(10) 



(11) 



We can see that equation ((TT|) is Fuchsian and it has four regular singular points Z{ 
over the Riemann sphere CP 1 . The infinity is an ordinary point for this equation. We 
assumed that k G (0,1), so s G (0, 00). Here it is important to notice that for real 
values of a, c, and d, and for all s G (0, 00) equation (JTTJl has four distinct regular 
singular points, i.e, the number of singular points does not depend on s. For further 
calculations we fix s = 1. 

Let us note here that transformation dTOb is a branched double covering of Riemann 
sphere CP 1 — ► 1\. Moreover, the branching points of this covering are precisely the 
four points where equation (fTTj) has singularities. 

Making the following change of the dependent variable 



m = w exp 



20 
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we can simplify ((TTJ to the standard reduced form 



w" = r{z)w, r(z) = lp'(z) + \ V {zf - q(z), (12) 



2 r v ' 4 

where r(z) can be written as 



r(z) = £ 



4 ' ' + * 



Z;')^ Z Zj 



(13) 



with coefficients 



«2 = = F — — + id, «3 = «4 = — , 

&l=&2 = d+i^F-|V 

5 j . c 7 



The differences of exponents A, = y/T+ 4fl^ at singular point Zj are the following 

Aj = A^ = VF - id, A 3 = A 4 = |. (14) 

3.2 Case L 3 = 

Let us assume that L3 = 0. Then, obviously d = 0, and, as we already mentioned, 
the special frame coincides with the principal axes frame. As we consider a symmetric 
body with A = B = 1, then additionally we have / = 0, and 

F= £_Z. _1 
2 16' C C 

Thanks to that equation (fl2)l has a simpler form and it can be transformed to a Riemann 
P equation. Instead of making a direct transformation in ((12)) it is more convenient to 
start from equation ©. Then, instead of transformation dlOt we make the following 
one 

f,-> z:=Ni{t,k) 2 , (15) 

and we obtain 

d 2 1/1 1 \ d 1-c/l 1 



, m + ~ r + ^" -n*w + — s— t m = 0. (16) 

d£ 2 2\z-l 2z/dt 8 V 2 " 1 z / 

For this Riemann P equation the difference of exponents at z = is 3/4, at z = 1 is 
1/2, and at z = 00 is 

Aoo = \V^7. (17) 

Let us notice that as in the case d / transformation dl~5b is a branched covering of 
Riemann sphere CP 1 — > Y^, and the branching points of this covering are precisely the 
three points where equation (fl6)l has singularities. 
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3.3 Case L3 = 0. Second particular solution. 

When L 3 = we have at our disposal another family of particular solutions. Under 
our assumption A — B — 1 and L = [1, 0, 0] T , the following manifold 

M := {(M, N) G C 6 I Mi = M 2 = N 3 = 0,N 2 + N| = 1} C M , 

is invariant with respect to the flow of system ©• Similarly as J\f , manifold M\ is 
diffeomorphic to T§^ C TS^. and it is a symplectic sub-manifold of Mo- The Euler- 
Poisson equations restricted to M\, read 



%u 3 



d d 

-N 2 , ^Ni = cM 3 N 2/ -gN 2 



-CM3N1. (18) 
We consider a family of phase curves k i—» of the above equations given by 

(19) 



1 9 

-cMf + Ni = e, 



n 2 + n| = i, 



where e = 2fc 2 - 1. For k e (0, 1) curves are non-degenerate elliptic curves. Varia- 
tional equations along phase curve have the following form 





mi 







(c-l)M 3 











0" 




mi 




ni2 




(1 - c)M 3 














1 




m 2 


d 


m 3 
















-1 







m 3 


df 


ni 










cN 2 





cM 3 







ni 




n 2 










-cNi 


-cM 3 










n 2 




"3 




-N 2 


Ni 
















"3 



and they have the following first integrals 

h = cM^m^ + ni, 
hi = Ninii + N 2 m 2 + M^n^, 
h 2 = Nifti + N 2 n 2 . 

The normal variational equations are given by 

_ Wl = (c - l)M 3 m 2/ 
^-m 2 = (1 - c)M 3 mi + n 3/ 



-N 2 nii + Nim 2 , 
= Nimi + N 2 m 2 + M 3 n 3 . 
Now, the reduction of the above system to the second order equation gives 



fi + (Ml - Ni)n = 0, 



n 



n 3 . 



(20) 



Let us notice that from equations dl8b and dl9b it follows that 

n -2c(e-Ni)(l-Nl). 

12 



Ni 2 



Thus, putting 

2 e 
Nt = -v+ -, 

c 3 

we obtain the following equation 

v 2 = Av 3 - g 2 v - g 3/ 
determining the Weierstrass function v(t) = p(t;g 2 ,g 3 ) with invariants 

g2 = lc 2 (e 2 + 3), g3 = l-c 3 e(e 2 -9). 



Hence, we can express N\, and M 2 (using (fl9)l ), in terms of the Weierstrass function 
p(t;g2,g3). The discriminant and the modular function of p(t;g 2r g 3 ) are following 

A = g 3 -27g 2 = c 6 (e 2 -l) 2 , 

Si _ (e 2 + 3) 3 



j(g2,g3) = j(e) 



g\-T7g 2 27(e*-l)2- 
Hence, we can rewrite equation d20l in the form of the Lame equation 

d 2 

= (ap{t;g 2 ,g 3 ) + P)n, (21) 

where 

a = 2C(2C + l), /3 = ^C(1-4C). 
It is important to notice here the physical restriction on parameter C, namely, we have 

ce (0,2). 



4 Proof of Theorem 

In our proof of Theorem |2] we try to be as close as possible to the proof of Ziglin. 
Namely, first we show that a necessary condition for integrability is L3 = (or L\ = 0, 
but this gives the already known integrable case of Lagrange). In fact this is the most 
difficult part of the proof. Then, we use the second family of particular solutions and 
we restrict the possible values of the principal moment of inertia. Finally, using the 
first solution, we limit all allowable values of C to those corresponding to the known 
integrable cases. 

We organise the three steps of the proof in the form of three lemmas. Only the first 
one is somewhat involved, the remaining two are very simple. 

The first step is to show that a necessary condition for integrability is d = 0, see 
formulae ©. 

Lemma 2. Let us assume that d 7^ 0. Then the identity component of the differential Galois 
group of equation ((T2)l is not Abelian. 
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Proof. In our proof we use of Lemma|8]and|9j see Appendix. If ((12)) is reducible then the 
identity component Q° of its differential Galois group is Abelian in two cases: when Q 
is a subgroup of diagonal group V or when it is a proper subgroup of triangular group 
7. 

First we show that Q (jL T>. Let us assume the opposite. Then there exist two 
exponential solutions of (|l2")l which have the following form 

4 

w x = Pj H(z - Zi ) e », P, € C[z], I = 1,2, 

i=l 

where e^i for 1 — 1,2 are exponents at singular point z„ i.e., 

eye ||(l + Ai),|(l-Ai)}. 

Here Aj for i = 1, ... ,4 are given by (fl~4|) . The product of these solutions v = W\Wj_ 
belongs to C(z) and it is a solution of the second symmetric power of ((12)) , i.e. equa- 
tion ((33)> with r given by ((l3)l . This equation has the same singular points as equa- 
tion (fl2)l . Exponents pq at singular points Zj, and at infinity p^i for the second sym- 
metric power of (fl2)l are given by 

p u e {1,1 ± Ai}, PoQ/ i e {-2,-1,0}, 1 = 1,2,3, 

where Aj for z = 1, 2, 3, 4 are given by (|14|| . If we write v = P/Q with P, Q G C[z] then 

K 

! = 1 

and n z = —pn G N for certain I. However, if d ^ 0, then p ; / is not a negative integer 
for z = 1,2,3,4 and I = 1,2,3. This implies that Q — 1. Hence, equation d33l has 
a polynomial solution v = P, and degP = — Poo,/ < 2. But v is a product of two 
exponential solutions of the form ||4jl, so we also have 

4 

v = PiP 2 n( z - Zi) ei ' m+ey G C [z]/ m > 1 G I 1 ' 2 )- C 22 ) 

i=l 

Consequently, + e z y is a non-negative integer for z = 1, 2, 3, 4. As for d ^ 0, we have 
2e z y Z for z = 1, 2, 3, 4 and I = 1,2, we deduce that in ((22)l z« 7^ Z. But e,-^ + 2 = 1/ 
for z = 1, 2, 3, 4. Thus, we have 

degv = degP = deg(PiP 2 ) + 4 > 4. 

We have a contradiction because we already showed that deg P < 2. 

It is also impossible that Q conjugates to T m for a certain meff because when d 7^ 
exponents for zj and z 2 are not rational. This implies also that Q is not finite. 

The last possibility that Q° is Abelian occurs when Q is conjugated with a subgroup 
of D\ We show that it is impossible. To this end, we apply the second case of the 
Kovacic algorithm, see Appendix. The auxiliary sets for singular points are following 

Ex = E 2 = {2}, E 3 = E 4 = {-1,2,5}, = {0,2,4}. 
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In the Cartesian product E = E M x Ej x • • • x E4 we look for such elements e for which 




is a non-negative integer. There are two such elements, namely 

e« = (4,2,2,-1,-1), = (2,2,2,-1,-1). 

We have d(e^) = 1 and d(e^) = 0. We have to check if there exists polynomial 
P = f\z + po which satisfies the following equation 

P'" + 39P" + (39 2 + 39' - 4r)P' + (d" + 309' + 9 3 - 4r9 - 2r')P = 0, 

where 

1 1 11 11 

9 = 1 . 

z-i z+i 2z- 1 2z+ 1 

Inserting P into the above equation we obtain the following system of linear equations 
for its coefficients 

-dp + Fp 1 =0, (1 -2F)p + 6dp 1 = 
Fp + rfpi=0, -6dp + (1 -2F)pi = 0. 

The above system for po an d pi has a non-zero solution if d 2 + F 2 = 0, but for a real d 
and F it is possible only when d = and F = 0. □ 

As the covering t z given by (flO)l does not change the identity component of 
the differential Galois group of the normal variational equations ©, from the above 
lemma it follows that if the Euler-Poisson equations are integrable, then 

For a symmetric body when A = B 7^ C, the above condition implies that either L\ = 0, 
and this corresponds to the integrable case of Lagrange, or L3 = 0. Hence, we have to 
investigate the last case. Notice that in this case the special frame is the principal axes 
frame so we have L = L = [1, 0, 0] T , a = 1/A = 1/B = b = 1 and c = 1/C ^ 1. At this 
point it is worth to observe that now the identity component of the differential Galois 
group of equation ((l2)l is Abelian for infinitely many values of c. Let us remind here 
that from the physical restriction it follows only that c G (1/2, 00). 

Proposition 1. If L3 = then the identity component of the differential Galois group of 
equation ((T2)l is Abelian in the following cases: 

c = 1 + 21(41 -1), c = 4 + 21(41 + 5), c = 2 + 21(41 + 3,) (23) 

11 79 4 

c=— +21(1 + 1), c =- + -l + 2l 2 , (24) 

where I is an integer. 
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Proof. When L3 = then the identity component of the differential Galois group of 
equation ((12)) is Abelian if and only if the differential Galois group of equation dl6t is 
Abelian (transformation from ((12)| to (fT6t is algebraic). Then applying Kimura Theo- 
rem |5] to equation ((l6)l we easily derive the above values of c for which the identity 
component of the differential Galois group of this equation is Abelian. □ 

Thus, applying the Morales-Ramis or Ziglin theory and using the first particular 
solution we cannot prove non-integrability of the Euler-Poison equations for all values 
of c listed in the above proposition. This is why, in the lemma below, we consider 
normal variational equations corresponding to the second particular solution. 

Lemma 3. Assume that C G (0,2) and C 7^ m/4, for m = 1,...,7, then for almost all 
eGl, the differential Galois group of equation ||2TJ is SL(2, C). 

Proof. We assume first that e ^ ±1. Then the discriminant of the elliptic curve associ- 
ated with p(t;g2,g3) does not vanish, and we can apply Lemma llOl see Appendix. We 
consider successively three cases from this lemma. 

For the Lame-Hermite case, we have a = n(n + 1) for n G Z. This implies that 
C = n/2, and hence, as C <E (0,2) we have C G {1/2,1,3/2}. 

For the Brioschi-Halphen-Crawford case, we have ot — n{n + 1) and ra = n + l/2G 
N. Thus, we have 

1 1 

C = -- + -m, m G N. 
4 2 

So, this case can occur only when C G {1/4,3/4,5/4,7/4}. 
In the Baldassarri case we notice that the mapping 

R\{-l,l}3e^j(e), 

is non-constant and continuous. Hence, by Dwork Proposition|2j for a fixed C this case 
can occur only for a finite number of values of e. □ 

As C = 1/2, C = 1 and C = 1/4 correspond to the Kovalevskaya, Euler and 
Goryachev-Chaplygin cases, respectively, we have to investigate cases 

CG {3/4,5/4,3/2,7/4}. 

To this end we return to equation (|T2j. As we show, for d = it can be transformed to 
the form dl6t and, moreover, this transformation does not change the identity compo- 
nent of its differential Galois group. We can prove the following. 

Lemma 4. For C G {3/4, 5/4, 3/2, 7/4} the differential Galois group of is SL(2, C). 

Proo/. For CG {3/4, 5/4, 3/2, 7/4} the respective values of the difference of exponents 
at infinity Aoo (see formula (|1Z|>) for equation ((16)) are following 

if /TT . /3 . [5 . fW] 

Now, a direct inspection of possibilities in the Kimura Theorem |5] shows that Riemann 
P equation (fl6)> with prescribed differences of exponents does not possess a Liouvillian 
solution, so its differential Galois group is SL(2, C). □ 

Now the proof of Theorem ||2jl is a simple consequence of the above three lemmas. 
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5 Proof of Theorem |3l 

In the proof of Theorem |3] we apply the Ziglin Lemma [l] and his idea of its application 
given in |52f| . 

The Euler-Poisson equations restricted to J\f possess a hyperbolic equilibrium at 
u = (0,0,0,1,0,0). The phase curve l~i corresponds to the solution of equations 
with k = 1. It contains two real components which are real phase curves corresponding 
to real solutions homoclinic to u. Their union is Re V\ and we denote its closure by CI. 

Lemma 5. For an arbitrary complex neighbourhood UcAfofCl there exists e > 0, such that 
for < 1 — k < e the fundamental group 7ii(V k ) of phase curve V k is generated by loops lying 
in U. 

Proof. The time parametrisation of is given by 

M 2 (t,k) := -2fccn(t,fc), 

Ni(t,k) := 2k 2 sn 2 (t,k) - 1, (25) 
N 3 (t,k) := 2fcsn(t,fc) dn(t,k), 

where sn(f,fc), cn(f,fc) and dn(t,k) denote the Jacobi elliptic functions of argument t 
and modulus k. Thus, particular solutions of (|lj 

<p{t,k) := (0,M 2 (U),0,Ni(t,fc),0,N 3 (f,fc)), (26) 
defined by d25b are single-valued, meromorphic, and double periodic with periods 

T x {k) = 2K(k) + 2iK'(k), T 2 (k) = 2K(k) - 2iK'(k), (27) 

where K(k) is the complete elliptic integral of the first kind with modulus k, K'(k) := 
K(k'), and k' := Vl — k 2 . In each period cell they have two simple poles at: 

h(k) = iK'(k), t 2 (k) = -iK'(k) mod (T 1 (k),T 2 (k)). 

Periods T\(k) and T 2 (k) given by <(27)> of solution ((26)l are primitive. Minimal real and 
imaginary periods are T(k) = 4K(k) and T'(k) = 4iK'(k). As abase point x(k) € 1\ we 
choose x(k) = (p(to(k),k) where to(k) = K(k). Let us notice that from ((25)) it follows 
that 

M 2 (t Q (k),k)=0, 

N l (t Q (k),k)=2k 2 -l, (28) 
N 3 (t (k),k)=2kk'. 

Now, let 

A k ,A' k : [0,1]-* T k , 

be the loops with base point x(k) corresponding to periods T(k) and T' (k), respectively. 
These loops cross at point 

x'(k) = <p(t (k) + T(k)/2,k) = <p(t (k) + T'(k)/2,k). 
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Figure 1: Parallelogram of period with marked loops 



As a results, we obtain four semi-loops with end points x(k) and x'(k). The funda- 
mental group TC\{T^,x{k)) of V% is generated by these semi-loops, see Figure [U Let us 
analyse what happens when k tends to 1. From d28l it follows that x(k) tends to u and 
from d25l we deduce that loop \ tends to £1. To see what happens with loop A[ when 
k tends to 1, let us put t = to(k) + ir in formulae (|25)l . We obtain 

Nl(a) =- 1 + 5nfe 
V 7 dn 2 (T,it') 

Thus, loop tends to point u ask tends to 1. □ 
Time parametrisation of these phase curves is given by 



2k 

M 3 (t,k) = — cn( cut, k), 
cu 

N!(t,k) = 2k 2 sn 2 (tut, k) - 1, 
N 2 (t,k) =2ksn(cut,k)dn(cut,k), 

where cu 2 = c. Thus, the second family of particular solutions of © 

<pi(U) := (0 f f M 3 (t f k) / N 1 (t f k) f N 2 (t f k) / 0) / 



(29) 
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defined by <(29|> contains solutions which are single-valued, meromorphic, and double 
periodic with periods 

Mk) n .K'(k) , m K(k) n .K'(k) 
T^k) =2-^- + 2i—^-, T 2 (k) = 2^^- - 2i— 

In each period cell they have two simple poles at: 

*!(*)=£££., = mod (T x (fc),T 2 (fc)). 

Let us notice that for the Euler-Poisson equations restricted to M\ point u = (0, 0, 0, 1, 0, 0) 
is also a hyperbolic equilibrium. The phase curve corresponds to the solution of 
equations (fl8)> with k = 1. As in the previous case, it contains two real components 
which are real phase curves corresponding to real solutions homoclinic to u. Their 
union is Re Vj and we denote its closure by Qj. Using the same arguments as in the 
proof of Lemma |5] we can show the following. 

Lemma 6. For an arbitrary complex neighbourhood U C Af\ of Cli there exists e > 0, such 
that for < 1 — k < e the fundamental group 7ti(V£) of phase curve 1\ is generated by loops 
lying in U. 

Now, to prove Theorem |3] let us notice that we showed that, except for the known 
integrable cases, the identity component of the differential Galois group of the normal 
variational equations corresponding to 1\ or is not Abelian for almost all values 
of k £ (0, 1). In fact, in Lemma 12 we proved that for d ^ the identity component 
of the differential Galois group of the normal variation equations corresponding to 1\ 
with k — 1/ \[2 is not Abelian. By Lemma it is not Abelian for almost all values of 
k G (0, 1). Then, in Lemma HI we proved that the identity component of the differential 
Galois group of the normal variation equations corresponding to is not Abelian for 
almost all values of k G (0, 1), except for C = m/4, m — 1, . . . , 7. Finally, for C = ra/4 
such that C ^ {1/4, 1/2, 1}, we showed that the identity component of the differential 
Galois group of the normal variation equations corresponding to with k = 1 / y/2 is 
not Abelian. Again, by Lemma|3it is not Abelian for almost all values of k € (0, 1). 

Both normal variational equations corresponding to 1\ and are Fuchsian. For 
a Fuchsian equation we know that if the identity component of its differential Galois 
group is not Abelian then its monodromy group does not possess a rational invariant, 
see Theorem 3.17 in [4]. 

Assume now that for C ^ {1/4,1/2,1} the Euler-Poisson equations possess an 
additional real meromorphic first integral defined in a real neighbourhood of flU D.\. 
Then we can extend this integral to a complex meromorphic one, defined in a certain 
complex neighbourhood U of CI U D.\. Then, by Lemma |5] and |6j we find such e > 
that the fundamental groups of F/ c and with < 1 — k < e, are generated by loops 
lying entirely in U. Then, from the Ziglin Lemma [TJ it follows that both monodromy 
groups of normal variational equations corresponding to 1\ and l~£ possess a rational 
invariant. However, above we showed that at least for one of them it is not true. A 
contradiction proves Theorem |3l 
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6 Remarks and Comments 



One important difference between the Ziglin and Morales-Ramis theory is related with 
the procedure of obtaining the normal variational equations. Assume that system (H} 
possesses certain number of known first integrals Hj, such that their differentials dff,, 
are linearly independent on F. Then dff, o 7T _1 for i = 1, . . . ,k are independent first 
integrals of ©.Their common level 

N p := {77 g F I dHj o n' 1 !] = p if p f <E C, z = 1, . . . , fc}, 

defines a ra-dimensional linear bundle over F, where m = n — k — 1. Using these 
integrals we can reduce the order of system (JSJ). Namely, we consider the reduced 
normal variational equations 

fj = 71*. (T(o)7T~* 77), vj e N p . (30) 

However after this reduction defined by Ziglin, instead of a linear, we have an affine 
bundle over F, equations d30l are generally not homogeneous ones, and the mon- 
odromy group is a subgroup of affine transformations of C m . Till now this construction 
has not been translated to the Morales-Ramis theory where we work with a system of 
homogeneous equations defined on No- To realise the importance of Ziglin reduc- 
tion let us notice that using only the reduced equation on No it is impossible to prove 
global non-integrability of the Goryachev-Chaplygin case. In [52] he gave such a non- 
integrability proof investigating the reduced normal variational equations on Np with 

The Morales-Ramis theory is coordinate independent, however, investigating a spe- 
cific problem, we always have to choose appropriate coordinates. The form of normal 
variational equations depends on local coordinates and this is why their choice is im- 
portant. It is especially evident when we investigate a problem connected with a rigid 
body. Equations of motion of the heavy top can be written in many different forms. 
As we mentioned in Remark [II the natural phase space for a rigid body with a fixed 
point is T*SO(3, R). There are no 'natural' coordinates on SO (3, R), and thus there are 
no 'natural' canonical coordinates. The most widely used are Androyer-Deprit canon- 
ical coordinates [8] or the Euler angles and conjugated momenta. In fact we checked 
which, from almost all known coordinates on T*SO(3, R), are most feasible for appli- 
cation of the Morales-Ramis theory. In our exposition we work with the Euler-Poisson 
equations. However, our choice of the body fixed frame is not conventional. Usually 
the principal axes frame is used. To see what is an advantage of our choice, let us no- 
tice that using the principal axes frame we can derive the normal variational equation 
in the form similar to ©, however, to put it in the form of an equation with ratio- 
nal coefficient we have to choose a transformation different than (ITOl , and as a result, 
we obtain, instead of equation ((12)) possessing four regular singularities, a much more 
complicated Fuchsian equation with seven singular points. Our choice of the body 
fixed frame appears e.g. in fiUl . 

Simplifications of the normal variational equations which occur when L3 = need 
an explanation. In fact, one can observe that although the Riemann surface F^ for the 
second particular solution is a torus with two points removed, see formulae <(29|> . the 
normal variational equation (|2l|) corresponding to it has the form of a Lame equation, 
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so it is defined over a torus with one point removed. The reason of what happened 
is symmetry. When L3 = the Euler-Poisson equations restricted to Mo are invariant 
with respect to an involutive symplectic diffeomorphism J\ : Mo 1— ► Mo defined by 

J 1 (M 1/ M 2/ M 3/ N 1/ N 2/ N 3 ) = (-Mi,M 2 ,-At3,Ni,-N2,N3). 
Let us denote 

M = {p G M I Ji{p) ^p}, M = M/J lt 

and let n : M 1— > M be the projection. In the natural way equations ((TJl induce Hamil- 
tonian equations on M with Hamiltonian function H = H o 7T _1 . Then, according 
to Ziglin, see Lemma on page 36 in IqO'], if system ((TJ is integrable, then the induced 
Hamiltonian system on M is also integrable. For the induced system we have a family 
of particular solutions (pi(t,k) = no cpi(t,k). The corresponding Riemann surfaces r£ 
are tori with one point removed. 

A simplification of the normal variational equation for solution (p(t,k) when L3 = 
and the fact that we can transform them to a Riemann P equation (for an appropriate 
choice of energy) is also related with symmetry. Namely, when L3 = 0, system (|TJ) 
restricted to Mo is also invariant with respect to an involutive symplectic diffeomor- 
phism J : Mo 1 — ► Mo defined by 

J (Mi, M2, M3, Ni, N2, N3) = (-M 1/ -M 2/ M3 / N l/ N2 f -N 3 ). 

For symmetry reduction of variational equations see Section 4.2 in [4]. 

Let us note that in Lemma 12 we claim that if d 7^ then the identity component of 
the differential Galois group of equation (fl2)l is not Abelian. Thus, it can be the whole 
group SL(2, C) or whole triangular subgroup T of SL(2, C). We do not know if the 
second case can occur. 

In the case of first particular solution we do not work with the elliptic curve 1\ but 
with the Riemann sphere (minus singular point) for which 1\ is a covering. The reason 
of this is that we have no tool similar to the Kovacic algorithm for a second order linear 
differential equation defined on an elliptic curve. However, in the case of the second 
particular solution we can work directly on elliptic curve because, in this case, the 
normal variational equation is the Lame equation for which the monodromy group is 
know. Of course, in this case we can also work on the Riemann sphere making well 
know transformation of the Lame equation to its algebraic form. 

In his proof of Theorem[2]and[3]Ziglin used the explicit time parametrisation of par- 
ticular solutions. First he showed that if the system is integrable then the monodromy 
of the normal variational equations along real periods of a particular solution must be 
equal to the identity. Then, using analytical tools he derived the necessary conditions 
for the integrability. 

In our exposition we use the explicit time parametrisation of particular solutions 
in the proof of Theorem |3l In fact, we use it only to show explicitly what happens 
with loops along real and imaginary periods when k tends to 1. However, one can 
deduce this information from the equations defining the elliptic curve. Thus, we can 
avoid using explicit time parametrisation at all. We keep it in the proof of Theorem |3] 
because, as we hope, it makes the exposition more transparent. 
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The physical restriction C G (0, 2) plays crucial role in our, as well as, in the Ziglin 
proof. Integral? le systems are really rare, hence it is an interesting question if the Euler- 
Poisson equations are integrable for values of parameters which do not satisfy this 
restriction. 

Considering the case L3 = and the first particular solution Ziglin showed that 
the necessary condition for integrability is c G N. In our Proposition [l] there are two 
families of c such that c G Q. The reason why they appears is that we fixed the energy 
for the first solution. 

In the proof of Lemma |3] we show that the Brioschi-Halphen-Crawford case is pos- 
sible only when C G {1/4,3/4,5/4,7/4}. For this values of C we can calculate the 
Brioschi determinant Q m (g2,g3, fi) defined by (|36)l . Calculations show that it vanishes 
identically only when C = 1/4, i.e., for the Goryachev-Chaplygin case. Thus, in fact, 
to prove that for C G {3/4,5/4,7/4} the Euler-Poisson equations are non-integrable, 
we can use the second solution. We use the first one because calculations are simpler. 
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8 Appendix 

8.1 Dependence on a parameter 

Let us consider a second order differential equation of the following form 

y" = r{z,e)y, ' = ^, (31) 

where r(z, e) is a rational function with respect to z and e, i.e., r G C(e) (z) = C(z, e). 
Here e plays the role of a parameter. For a fixed value of e we denote by G°(£) the 
identity component of the differential Galois group of equation (|3l) . Let U C C denote 
an open not empty connected set with compact closure. We show the following. 

Lemma 7. Assume that: 
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1. Equation (|3~T]l is Fuchsian. 

2. For e £ IT, equation d3Tb possesses N singular points (N does not depend on e)for which 
exponents do not depend on e. 

3. For eo £ LL {?°(£o) & noi solvable (is not Abelian). 

Then, except finitely many values ofe^U, Q°(s) is not solvable (is not Abelian). 

Proof. From the Kovacic algorithm in the form given in fl3l: l38|l we know that, under 
our assumption, if G°(e) is solvable, then there exists a polynomial P (whose degree 
does not depend on e) which is a solution of a linear differential equation L(y) = 
with coefficients in C(z, e). The order of L(y) = does not depend on e. We have a 
finite number of choices for the degree of P and a finite number of choices of L(y) = 0. 
Finding a polynomial solution of linear equation L(y) = reduces to finding a non- 
trivial solution of a homogeneous linear system with coefficients in C(e). But the last 
problem reduces to finding common zeros of a finite number of polynomials. We know 
that not all of these polynomials vanish identically (otherwise Q q {eq) is solvable). Thus, 
there is at most a finite number of values of e for which they vanish simultaneously. 
Finally, let us notice that set 

{e £ U | g°(e) is Abelian}, 

is a subset of 

{e £ U\ g°(c) is solvable}. 

□ 



8.2 Second order differential equations with rational coefficients 

Let us consider a second order differential equation of the following form 



y" = ry, 



r e C(z), 



d_ 
dz 



(32) 



For this equation its differential Galois group Q is an algebraic subgroup of SL(2, C). 
The following lemma describes all possible types of Q and relates these types to forms 
of solution of 63, see OB- 

Lemma 8. Let Q be the differential Galois group of equation <(32|> . Then one of four cases can 
occur. 

1. Q is reducible (it is conjugated to a subgroup of triangular group) ; in this case equation 
d32b has an exponential solution of the form y = exp / tu, where cu £ C(z), 



2. Q is conjugated with a subgroup of 



Dt 



c 

c" 1 



c e C* u 



c 
r 1 



c £ 



in this case equation ((32)l has a solution of the form y = exp / tu, where cu is algebraic 
over C(z) of degree 2, 



23 



3. Q is primitive and finite; in this case all solutions of equation (|32)l are algebraic, 

4. Q = SL(2, C) and equation d32b has no Liouvillian solution. 

We need a more precise characterisation of case 1 in the above lemma. It is given 
by the following lemma, see Lemma 4.2 in Hzi . 

Lemma 9. Let Q be the differential Galois group of equation d32b and assume that Q is re- 
ducible. Then either 

1. equation d32b has a unique solution y such that y'/y G C(z), and Q is conjugate to a 
subgroup of the triangular group 



T 



a b 
a' 1 



a,b G C,a ^ 0} . 



Moreover, Q is a proper subgroup of T if and only if there exists m G N such that 
yin ^ C(z)_ In this case Q is conjugate to 



T 



a b 
a' 1 



a, be C,a m = 1 



where m is the smallest positive integer such that y m G C(z), or 

2. equation d32b has two linearly independent solutions y\ and yi such that y'Jyi G C(z), 
then Q is conjugate to a subgroup of 



V 



a 
a' 1 



In this case, y\yi G C(z). furthermore, Q is conjugate to a proper subgroup ofV if and 
only ifyf G C(z) for some m G N. In this case Q is a cyclic group of order m where m 
is the smallest positive integer such that yf G C(z). 

In case 2 of the above lemma we know that v = yij/2 G C(z). Differentiating v three 
times, and using the fact that y ( satisfies equation <(32|> . we obtain 



./// 



2r'v + 4rv'. 



(33) 



The above equation is called the second symmetric power of equation d32b . For applica- 
tions of symmetric powers of differential operators to study the existence of Liouvillian 
solutions and differential Galois group see e.g. 1 42 ; S43l l46f| - 

To decide if case 2 from Lemma |8] occurs we can apply the Kovacic algorithm. Here 
we present its part devoted to this case and adopted to a Fuchsian equation. At first 
we introduce notation. We write r(z) G C(z) in the form 



r(z) 



sjz) 
tiz) 



s(z), f(z) G C[z], 



where s(z) and t(z) are relatively prime polynomials and t(z) is monic. The roots of 
t(z) are poles of r(z). We denote I' := {c G C | t(c) = 0} and I := I' U {oo}. The 
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order ord(c) of c G L' is equal to the multiplicity of c as a root of t(z), the order of 
infinity is defined by 

ord(oo) := max(0,4 + degs — degf). 

Because we assume that equation (|32)l is Fuchsian, we have ord(c) < 2 for c £ I. 
For each c G L' we have the following expansion 

and we define A c = \/\ + ^a c . For infinity we have 

r(z) = ^ + o(iV 



z 2 \z 3 

and we define Aoo = ^/l + ^a^. 

The algorithm consists of three steps. 
Step I. For c G L 1 such that ord(c) = 1 we define E c = {4}; if ord(c) = 2 

E c := {2,2(1 + A c ),2(l-A c )}nZ. 

If ord(oo) < 2 we put Eqo = {0,2,4}; if ord(oo) = 2 we define 

:= {2,2(1 + A oc ),2(l - A^)} H Z. 

Step II. For each e in the Cartesian product 

E := Eqo x J~J E c , 

eel' 



we compute 



d{e) := ^ ( - e t 



eel' 

We select those elements e G E for which rf(e) is a non-negative integer. If there are no 
such elements Case 2 from Lemma |8]cannot occur and the algorithm stops here. 
Step III. For each element e G E such that d(e) = n G No we define 

ceL' 

and we search for a monic polynomial P = P(z) of degree n satisfying the following 
equation 

P'" + 39P" + (39 2 + 39' - 4r)P ! + (d" + 399 1 + 9 3 - 4r9 - 2r')P = 0. 

If such polynomial exists, then equation <(32|> possesses a solution of the form w = 
exp / cu, where 

11 P 1 
w 2 + il)cu + + -ip 2 - r = 0, if) = 9+— . 

If we do not find such polynomial, then case 2 in Lemma El cannot occur. 
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8.3 Riemann P equation 



The Riemann P equation [48] is the most general second order differential equation 
with three regular singularities. If we place, using homography, these singularities at 
z = 0, 1, oo, then it has the form 

d 2 £ , fl-a-a' , l-y-y'\ d£ 

^ 04) 



dz 2 



+ 



+ 



aa 
+ -*- + 



z z — 1 

yy 7 1313' - aa! - yy' 



l) 2 



+ 



z z-1 



5 = 0, 



where (a, a'), (y, y') and (fl, (3') are the exponents at singular points. Exponents satisfy 
the Fuchs relation 

a + ot + y + y> + $ + &' = \. 
We denote differences of exponents by 

A = a — a' , -v = y — y' , \jl — (3 — (3' . 

For equation d34b the necessary and sufficient conditions for solvability of the identity 
component of its differential Galois group are given by the following theorem due to 
Kimura Ol. see also R 

Theorem 5 (Kimura). The identity component of the differential Galois group of equation 04b 
is solvable if and only if 

A: at least one of four numbers A + \x + y, —A + fx + y, A — \i + y, A + ju — v is an odd 
integer, or 

B: the numbers A or — A and \x or —\xand y or —y belong (in an arbitrary order) to some of 
the following fifteen families 



1 


1/2 + 1 


1/2 + s 


arbitrary complex number 




2 


1/2 + 1 


1/3 + s 


1/3 + q 




3 


2/3 + 1 


1/3 + s 


1/3 + q 


l + s + q even 


4 


1/2 + 1 


1/3 + s 


1/4 + q 




5 


2/3 + 1 


1/4 + s 


1/4 + q 


l + s + q even 


6 


1/2 + 1 


1/3 + s 


1/5 + q 




7 


2/5 + 1 


1/3 + s 


1/3 + q 


l + s + q even 


8 


2/3 + 1 


1/5 + s 


1/5 + q 


l + s + q even 


9 


1/2 + 1 


2/5 + s 


1/5 + q 


l + s + q even 


10 


3/5 + 1 


1/3 + s 


1/5 + q 


l + s + q even 


11 


2/5 + 1 


2/5 + s 


2/5 + q 


l + s + q even 


12 


2/3 + 1 


1/3 + s 


1/5 + q 


I + s + q even 


13 


4/5 + 1 


1/5 + s 


1/5 + q 


l + s + q even 


14 


1/2 + 1 


2/5 + s 


1/3 + q 


l + s + q even 


15 


3/5 + 1 


2/5 + s 


1/3 + q 


l + s + q even 



Here I, s and q are integers. 

The solvability conditions are sufficient for our purposes because if Q° is not solv- 
able, then obviously it is not Abelian. 
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8.4 Lame equation 

The Weierstrass form of the Lame equation is following 

d 2 y 



d£ 2 



(35) 



where a and j3 are, in general, complex parameters and p(t; gi, g3) is the elliptic Weier- 
strass function with invariants gi, g3- In other words, p{t;g2,g?>) is a solution of the 
differential equation 

v 1 = f(v), f(v) = 4v 3 - g 2 v - g 3 . 
It is assumed that equation f(v) = has three different roots, so 

A = g 3 2 -27g 2 3 ^0. 

We recall that the modular function j{g2,g?>) associated with the elliptic curve u 2 = 
4v 3 — g2V — g3 is defined as follows 



si 



g\ 



Classically the Lame equation is written with parameter n instead of a related by the 
formula a = n(n + 1). We see that the Lame equation depends on four parameters 
(n, fi/g2/g3)- The following lemma lists all the cases in which the identity component 
of the differential Galois group of Lame equation d35b is Abelian, see H38L Sec. 2.8.4]. 

Lemma 10. The identity component of the differential Galois group of Lame equation d35t is 
Abelian only in the following cases: 

1. the Lame-Hermite case when n G Z and three other parameters are arbitrary. 

2. the Brioschi-Halphen-Crowford case for which m := n + j G N, and remaining param- 
eters (g2,g3, /3) satisfy an algebraic equation 

Qm (g2,g3,P) = 0/ 

see below. 

3. the Baldassarri case 2n g" Z, and 

n±qeZ for some q £ {1/4,1/6,1/10,3/10}, 
with additional algebraic restrictions on (g2,g3, &)• 

Polynomial Qm{g2>g3> fi) which appears in the Brioschi-Halphen-Crowford case, 
called the Brioschi determinant, is defined as follows 










m — 1 











12,1 




2(m-2) 








93,1 


<?3,2 




3(m-3) 








^4,2 


<?4,3 




4(m - 4) 






















tfm—2,m—3 


C\m—\,m—2 













tfm,m—2 





j8 m- 1 

tfm,m—l 



(36) 
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where 

m+u = ^( 2m - 0( m - 0/ ii+2,i = ^( 2m - 0( 2m - * - !)■ 

Algebraic restrictions on {g2,g3, /3) in the Baldassarri case are involved. Instead of 
them we use the following proposition which follows from one unpublished result of 
B. Dwork, see B. 

Proposition 2. The Baldassarri case for equation <(35|> occurs only for a finite number of values 
fj(g2,g3)- 
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